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ABSTRACT 

As  a  part  of  the  Comprehensive  Test  Ban  Treaty  Research  and  Development  Program, 
we  have  started  to  evaluate  regional  seismic  discriminants  for  western  China.  Our  primary 
objective  is  to  appraise  the  current  seismic  discrimination  capability  for  central  and  south- 
central  Asia.  We  have  been  testing  several  discriminants  using  short-period  and  broadband 
waveforms  recorded  at  the  Chinese  Digital  Seismic  Network  station  WMQ.  Our  preliminary 
results  cover  four  classes  of  discriminants:  (1)  short-period  ratios,  (2)  short-period  spectrd 
ratios,  (3)  long-period  ratios,  and  (4)  long-period  to  short-period  ratios.  For  each  discriminant 
tested,  we  have  measured  waveforms  from  over  140  earthquakes  to  find  ratio-distance  rela¬ 
tionships,  and  we  have  used  these  relationships  to  correct  our  final  ratio-m^  discrimination 
plots  for  distance  effects. 

To  date,  we  have  evaluated  earthquakes  that  range  in  magnitude  (m^)  between  3.8  and 
6.2  and  occurred  up  to  1200  km  from  WMQ.  We  have  evaluated  about  15  underground 
nuclear  explosions  from  the  Kazakh  test  site  (4.5  <6.1  and  A =960  km)  and  a  single 

nuclear  explosion  from  the  Lop  Nor  test  site  =4.6  and  A =260  km).  For  this  range  of 
magnitudes  and  distances,  we  find  that  most  of  the  discriminants  outlined  above  (with  the 
exception  of  the  spectral  discriminant)  separate  the  explosion  and  earthquake  populations. 
We  find  that  the  absence  of  a  Love  wave  accompanied  by  the  presence  of  a  Rayleigh  wave  is 
a  strong  indicator  of  an  explosion.  Such  an  occurrence  is  apparently  rare  for  central  Asian 
earthquakes.  We  plan  to  continue  this  study  by,  processing  other  regional  data,  developing 
regional  rtif,,  M^,  and  coda  magnitude  scales,  and  testing  the  ntf,  /M^  and  coda  discriminants. 
We  will  be  applying  multivariate  techniques  to  combine  discriminants  and  assess  discrimina¬ 
tion  capabilities. 
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OBJECTIVE 

As  a  part  of  the  Comprehensive  Test  Ban  Treaty  (CTBT)  Research  and  Development 
Program,  we  have  started  to  evaluate  regional  seismic  discriminants  for  western  China.  Our 
primary  objective  is  to  appraise  the  current  seismic  discrimination  capability  for  central  and 
south-central  Asia.  We  have  been  testing  several  known  discriminants  (cf.  Pomeroy  et  al., 
1982;  Taylor  et  al,  1989)  using  waveforms  recorded  with  the  Chinese  Digital  Seismic  Net¬ 
work  (CDSN).  Here  we  present  results  for  station  WMQ  in  northwest  China.  These  prelim¬ 
inary  results  cover  four  classes  of  discriminants:  (1)  short-period  ratios  (such  as  Pg  ILg),  (2) 
short-period  spectral  ratios  (such  as  Lg  {1-2 Hz)  /  Lg  {6-^ Hz)),  (3)  long-period  ratios  (Love 
wave  (L)  /  Rayleigh  wave  (R),  and  (4)  long-period  to  short-period  ratios  (such  as  L  /P^).  We 
compute  distance  corrections  by  linear  regression  of  earthquake  ratios  on  the  logarithm  of  dis¬ 
tance  and  apply  the  corrections  to  both  earthquake  and  explosion  ratios.  To  date,  we  have 
analyzed  about  140  earthquakes  that  range  in  magnitude  (w^)  between  3.8  and  6.2  and  range 
in  epicentral  distance  between  about  100  and  1200  km.  We  have  evaluated  about  15  under¬ 
ground  nuclear  explosions  from  the  Kazakh  test  site  (KTS)  and  a  single  nuclear  explosion 
from  the  Lop  Nor  test  site.  The  explosions  range  in  magnitude  between  4.6  and  6.1  and 
range  in  epicentral  distance  between  about  260  and  980  km.  For  this  range  of  magnitudes  and 
distances,  we  find  that  most  of  the  discriminants  outlined  above  (with  the  exception  of  the 
spectral  discriminant)  separate  the  explosion  and  earthquake  populations. 

RESEARCH  ACCOMPLISHED 

Waveforms  obtained.  Our  data  source  for  this  project  has  been  the  IRIS  Data  Manage¬ 
ment  Center  (DMC).  Using  NEIC/ISC  event  catalogs  maintained  at  the  IRIS  DMC,  we  have 
requested  and  obtained  short-period  (SH,  40  sps)  and  broadband  (BH,  20  sps)  waveforms 
from  about  140  earthquakes  and  16  nuclear  explosions  recorded  at  station  WMQ.  The  loca¬ 
tions  of  these  events,  dong  with  station  locations  and  some  geographic  features,  are  shown  in 
Figure  1.  We  obtained  and  processed  data  from  events  shown  in  Ae  box  surrounding  station 
WMQ.  Waveforms  of  a  single  explosion  from  Lop  Nor  are  available  through  IRIS.  This 
explosion  is  about  260  km  from  WMQ  and  has  a  reported  of  4.6.  About  15  nuclear  explo¬ 
sions  from  KTS  have  been  recorded  at  WMQ.  The  KTS  explosions  are  between  about  950 
and  980  km  from  WMQ.  Magnitudes  of  these  explosions  range  between  4.6  and  6.1.  Most 
earthquakes  composing  our  data  set  occurred  to  the  south,  southwest,  and  west  of  WMQ. 
They  range  in  magnitude  between  3.8  and  6.2.  Maximum  epicentral  distance  for  the  earth¬ 
quakes  is  about  1200  km. 

Waveforms  and  Processing  -  Short  Periods.  To  evaluate  short-period  discriminants,  we 
measured  bandpass  filtered  velocity  records  for  P„,  Pg,  S„,  and  Lg  rms  amphtudes.  Prior  to 
measuring  phase  amplitudes,  we  handpicked  the  clearest  arrivals  for  each  phase  on  unfiltered, 
three-component  records,  and  measured  velocities  for  each  phase.  Using  the  velocities  of 
each  phase,  we  defined  group  velocity  windows  centered  on  each  phase,  and  then  automati¬ 
cally  measured  rms  amplitudes  within  the  windows.  The  centers  of  the  frequency  bands 
evaluated  range  between  1  and  9  Hz. 

An  earthquake  from  January,  1990  has  nearly  the  same  epicenter  and  as  the  Lop  Nor 
explosion  (Figure  1).  As  an  example  of  our  data  processing  procedure,  we  compare  short- 
period  records  of  this  explosion  and  earthquake  in  Figure  2.  Each  event  shows  the  unfiltered 
SHZ  velocity  record  and  two  bandpass  filtered  waveforms.  Also  shown  are  the  P  and  Lg 
measurement  windows.  P„  and  S„  were  not  measured  on  these  records  because  the  event- 
station  distance  is  too  short  for  a  clear  head  wave  arrival  from  the  crust-mantle  boundary. 
Note  that  for  the  1-2  Hz  band  the  P  and  L  amplitudes  of  the  explosion  are  about  the  same. 
However,  for  the  4-8  Hz  band,  the  P  amputude  is  much  greater  than  the  Lg  amplitude.  For 
the  earthquake,  the  Lg  amplitude  is  always  greater  than  the  Pg  amplitude. 

Waveforms  and  Processing  -  Long  Periods.  Figure  3  shows  Rayleigh  waves  and  Love 
waves  from  the  same  explosion  and  earthquake  discussed  above.  Here,  we  show  broadband 
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displacement  records  filtered  between  8  and  16  seconds.  By  visually  inspecting  many  long- 
period  records,  we  selected  L  and  R  group  velocities  that  allowed  us  to  automatically  pick 
peak-to-peak  maximum  amplitudes  for  each  surface  wave  phase.  We  also  processed  surface 
waves  for  6-12  s,  12-24  s,  and  16-32  s  periods.  L  and  R  from  the  earthquake  have  about  the 
same  maximum  amplitudes,  while  R  for  the  explosions  is  about  four  times  greater  than  L. 
Indeed,  L  from  the  explosion  source  is  such  a  weak,  poorly  developed  phase,  that  it  often 
barely  exceeds  noise  levels.  Although  these  events  have  the  same  ,  the  amplitudes  of  the 
earthquake  surface  waves  are  more  than  a  factor  of  10  greater  the  the  amplitudes  of  the  explo¬ 
sion  surface  waves. 

Short-Period  Discriminants.  Many  short-period  regional  discriminants  that  exploit  a.  P  I S 
ratio  have  been  discussed  in  the  literature  (cf.  Pomeroy  et  al.,  1982).  Figure  4  shows  two  of 
our  examples  from  China.  For  each  discriminant,  we  show  the  ratio-distance  trend  obtained 
from  earthquake  data,  and  we  present  a  ratio-m^  discrimination  plot  where  all  been  been 
corrected  for  distance.  All  ratios  shown  have  signal  that  is  at  least  2.5  times  the  measured 
noise  level.  The  1-2  Hz  PgILg  ratio  does  not  separate  the  explosion  and  earthquake  popula¬ 
tions  nearly  as  well  as  the  4-8  Hz  ratio.  We  find  that  for  frequency  bands  above  3  Hz,  the 
populations  can  be  separated  for  PgILg,  P^ILg,  and  Pn  ^ ratios.  However,  considering 
signal-to-noise  ratios  of  these  four  body  wave  pnases,  we  find  that  the  P„  I  Lg  ratio  is  most 
useful. 

Short-Period  Spectral  Discriminants.  Spectral  discriminants  have  been  successfully 
applied  to  relatively  small  events  from  the  western  United  States.  Taylor  et  al.  (1988)  used  the 
Lg{\—2Hz)  /  Lg(6-SHz)  ratio,  Bennett  and  Murphy  (1986)  used  the  (0.5-1  Hz)  / 
Lg{2-AHz)  ratio,  and  Hartse  et  al.  (1995)  used  the  coda  (0.5-1  Hz)/ (2-4 Hz)  ratio  to 
separate  nuclear  explosions  from  earthquakes.  However,  Figure  5  shows  that  the  spectral 
discriminant  for  western  China  may  not  be  as  successful.  The  Lg  (1-2 Hz)  /  Lg  (6— 8 Hz)  ratio 
fails  to  separate  the  events.  We  also  have  tried  Lg  (0.5-1  Hz)  /  Lg  (2-4 Hz)  ratios  and  vari¬ 
ous  Pg  /  Lg  spectral  ratios  with  results  similar  to  those  in  Figure  5.  The  western  U.S.  studies 
generally  examined  small  explosions  (m^  <4.8)  that  were  detonated  above  the  water  table. 
Most  explosions  from  Asia  are  thought  to  be  detonated  below  the  water  table  {Matzko,  1994). 
As  discussed  by  Taylor  and  Denny  (1991),  the  amount  of  gas-filled  porosity  near  the  explo¬ 
sion  appears  to  control  the  performance  of  the  spectral-ratio  discriminant. 

Long-Period  Discriminants.  We  tested  long-period  L  /  R  discriminants  in  four  frequency 
bauids,  and  found  that  the  8-16  s  period  shows  the  best  separation  and  least  scatter  (Figure  6). 
The  L  IR  discriminant  was  mentioned  by  Pomeroy  et  al.  (1982),  but  has  not  been  widely 
applied.  The  problem  is  that  the  weak  L  of  the  explosions  often  fail  to  exceed  noise  levels. 
The  Kazakh  explosions  had  to  be  large  (m^  >  5.5)  to  generate  L  that  was  above  the  noise  at 
WMQ,  over  950  km  from  the  station.  However,  for  smaller  events,  the  absence  of  L  com¬ 
bined  with  the  presence  of  R  can  be  used  as  an  explosion  indicator.  We  measured  L  and  R  of 
over  100  earthquakes  and  found  only  four  events  (<4%)  where  L  was  below  the  noise  while 
R  was  above  ^e  noise.  At  the  same  time,  of  the  15  Kazakh  explosions  we  processed,  we 
found  four  (>  25  %)  with  L  below  the  noise  and  R  above  the  noise. 

Short-PeriodI Long-Period  Discriminants.  We  tested  PgIL  and  PgIR  ratios.  These 
discriminants  are  analogous  to  the  discriminant  that  has  been  successfully  applied  in 

the  western  U.S  {Taylor  et  al.,  1989).  Because  we  have  not  yet  developed  regional  and 
scales  for  station  WMQ,  we  have  simply  tested  Pg  /  surface  wave  ratios.  We  measured 
peak-to-peak  displacement  amplitudes  from  broadband  records  to  test  this  discriminant.  Figure 
7  shows  the  Pg  (2-4 Hz)  /  L  (8— 16r)  and  the  P.  (2— 4 Hz)  /  R  (8-16^)  ratios  versus  m^.  We 
also  tested  Pg  (1-2 Hz),  but  found  that  Pg  (2-4 Hz)  produces  better  separation.  L  performs 
better  than  R ,  but,  as  in  the  case  of  the  long-period  discriminants,  L  for  smaller  explosions  is 
difficult  to  measure  above  the  noise. 
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RECOMMENDATIONS  AND  FUTURE  PLANS 

We  have  started  a  comprehensive  seismic  discrimination  study  of  western  China.  Using 
waveforms  from  over  100  earthquakes  and  16  nuclear  explosions  recorded  at  CDSN  station 
WMQ,  we  have  tested  short-period  phase  ratios,  short-period  spectral  ratios,  long-period 
ratios,  and  short-period/long-period  ratios  as  event  discriminants.  We  have  estimated  and 
applied  distance  corrections  for  every  discriminant  tested.  We  find  that  the  absence  of  a  Love 
wave  accompanied  by  the  presence  of  a  Rayleigh  wave  is  a  strong  indicator  of  an  explosion. 
This  appears  to  be  a  rare  occurrence  for  central  Asian  earthquakes.  With  the  exceptions  of 
the  short-period  spectral  ratio,  we  have  been  able  to  separate  earthquakes  from  underground 
nuclear  explosions  for  events  down  to  =4.5  at  distances  of  about  1000  km.  However,  we 
have  not  yet  processed  any  chemical  explosion  data.  We  do  not  know  what  our  discrimination 
capabilities  will  be  for  smaller  explosions. 

Our  study  has  progressed  to  where  the  performance  of  some  discriminants  can  be  com¬ 
pared  between  the  western  China  region  and  the  western  U.S.  In  Figure  8  we  compare  event 
separation  for  these  two  regions  using  the  Pg  I  Lg  ratio.  The  figure  shows  average  ratios  and 
one  standard  deviation  for  earthquakes  and  nuclear  explosions  plotted  against  frequency.  Com¬ 
pared  to  the  U.S.,  in  western  China,  separation  is  greater  and  occurs  at  lower  frequencies.  As 
more  events  are  processed  from  Asia,  the  separation  may  decrease.  However,  the  good  separa¬ 
tion  at  about  3  Hz  implies  that  this  regional  discriminant  can  be  applied  at  greater  distances  in 
central  Asia  than  in  the  western  U.S.  Additionally,  the  scatter  of  the  Pg  ILg  ratio  appears  to 
less  for  western  China  than  for  the  western  U.S.  at  higher  frequencies. 

Furture  work  will  involve  expanding  outward  from  western  China.  We  have  obtained, 
but  have  not  yet  processed,  waveforms  recorded  at  stations  AAK  and  TLY  (Figure  1),  and  we 
will  obtain  and  process  other  regional  data.  We  will  develop  regional  ,  and  coda  mag¬ 

nitude  scales,  and  test  the  ttij,  /M^  and  coda  discriminants.  Because  we  have  tested  several 
different  discriminants,  we  have  been  preparing  multivariate  techniques  that  combine  discrim¬ 
inants  and  discrimination  capabilities.  We  will  analyze  smaller  magnitude  events  by  acquiring 
Chinese  provincial  catalogs  and  obtaining  these  events  from  the  IRIS  DMC. 
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Figure  1.  Event  location  map  of  central  Asia  for  the  years  1989-1994.  Large  symbols  inside 
study  area  box  denote  events  for  which  we  have  obtained  and  processed  waveforms. 
Asterisks  denote  nuclear  explosions  and  circles  denote  earthquakes.  Small  dots  outside 
study  area  represent  earthquakes  that  we  have  not  investigated.  For  scale,  it  is  about  260 
km  from  Lop  Nor  to  station  WMQ  and  about  970  km  from  the  Kazakh  test  site  to  sta¬ 
tion  WMQ.  Locations  are  from  the  National  Earthquake  Information  Center  catalogs 
maintained  at  the  IRIS  Data  Management  Center. 
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Lop  Nor  Explosion  -  mb  =  4.6 
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Lop  Nor  Earthquake  -  mb  =  4.6 
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Figure  2.  Short-period,  vertical  component  (SHZ)  waveforms  recorded  at  station  WMQ  for 
the  September  29,  1988  Lop  Nor  nuclear  explosion  {top)  and  a  nearby  earthquake  from 
1990  {bottom)  (see  Figure  1  for  locations).  Unfiltered,  bandpass  filtered  from  1-2  Hz,  and 
bandpass  filtered  from  4-8  Hz  waveforms  are  displayed  in  velocity.  Pg  and  Lg  measure¬ 
ment  windows  are  marked  by  Pgl,  Pg2,  Lgl,  and  Lg2.  For  botii  events,  /n^  =4.6  and 
event-station  distance  is  about  260  km.  Note  the  strong  Pg  phase  of  the  explosions  and 
that  the  P-  /  L„  ratio  of  the  explosion  increases  as  passbands  are  increased  to  above  4 
Hz. 
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Lop  Nor  Explosion  -  Surface  Waves 


Lop  Nor  Earthquake  -  Surface  Wave 


Figure  3.  Lxjng-period  (8-16  s)  surfeice  waves  recorded  with  broadband  seismometers  at  sta¬ 
tion  WMQ.  These  are  the  same  events  shown  in  Figure  2,  but  the  long  periods  are 
displayed  in  displacement.  Top  shows  the  explosion  Rayleigh  wave  between  R1  and  R2 
markers  and  the  Love  wave  between  LI  and  L2  markers.  Bottom  shows  the  Rayleigh  and 
Love  wave  of  the  earthquake.  Note  that  the  Love  wave  for  the  explosion  is  extremely 
weak  relative  to  the  Rayleigh  wave,  while  the  Love  and  Rayleigh  wave  of  the  earthquake 
have  comparable  amplitudes.  Also,  despite  the  similar  magnitudes  (=4.6),  the  surface 
waves  of  the  explosion  are  much  weaker  than  those  of  the  earthquake. 
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LoglO  Pg  4-8  Hz  /  Lg  4-8  Hz  LoglO  Pg  1-2  Hz  /  Lg  1-2  Hz 
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Figure  4.  Short-period  discrimination  plots  for  station  WMQ.  Top  shows 
Pg  (1-2 Hz)/ (1-2 Hz)  results,  and  bottom  shows  Pg  (4-8 Hz) / Lg  (4-8 Hz)  results. 
Ratio-distance  trend  (slope  and  intercept)  is  obtained  from  the  earthquake  measurements, 
and  trend  is  removed  in  the  ratio-magnitude  discrimination  plots  (on  right).  Waveform 
rms  amplitudes  are  measured  from  short-period,  vertical  component  (SHZ)  velocity 
records.  Note  increased  separation  between  earthquakes  and  explosions  as  frequency 
increases. 
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Uncorrected  Lg  /  Lg  v  Diet  -  WMQ  -  0HZ  Corrected  Lg  /  Lg  v  Mag  -  TFMQ  -  $HZ 


Distance  -  km 


Figure  5.  Short-period  spectral  discrimination  plots  for  station  WMQ.  The  Lg  (\-2)/Lg  (6-8) 
spectral  discriminant  was  successfully  applied  to  events  from  the  western  U.S.  by  Taylor 
et  al.  (1988),  but,  for  south-central  Asia,  the  Lg  (l-2)/L^  (6-8)  spectral  ratio  apparently 
does  not  perform  well.  Distance  corrections  have  been  applied  to  the  ratio-magnitude  plot 
as  described  in  Figure  4.  Waveform  rms  amplitudes  are  measured  from  short-period, 
vertical-component  (SHZ)  velocity  records. 


Uncorrecled  L  /  R  v  Dist  -  WMQ  -  BHL  Corrected  L  /  R  v  Mag  -  WMQ  -  BHL 


Distance  -  km 


Figure  6.  Long-period  (8-16  s)  discrimination  plots  for  station  WMQ.  Rayleigh  wave  {R)  is 
measured  from  vertical  component,  and  Love  wave  (L )  is  measured  from  tangential  com¬ 
ponent  of  broadband  (BH)  displacement  records.  Distance  corrections  have  been  applied 
to  the  ratio-magnitude  plot  as  described  under  Figure  4.  When  Love-wave  energy 
exceeds  signal-to-noise  levels,  the  L  /  R  ratio  is  a  viable  discriminant. 
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LoglO  L  8-16  8  /  Pg  2-4  Hz  LoglO  R  8-16  s  /  Pg  2-4  Hz 


Uncorrected  R  /  Pg  v  Dist  -  WMQ  -  BLBZ  Corrected  R  /  Pg  v  Mag  -  WMQ  -  BLBZ 
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Figure  7.  Long-period  /  short-period  discrimination  plots  for  station  WMQ.  Top  shows 
R  (8-16s)/Pg  (2-4  Hz)  results,  and  bottom  shows  L  (8-16s)/Pg  (2-4  Hz)  results.  Dis¬ 
tance  corrections  have  been  applied  to  the  ratio-magnitude  plot  as  described  under  Figure 
4.  Waveform  peak-to-peak  amplitudes  are  measured  from  broadband  (BH)  displacement 
records. 
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Figure  8.  Short-period  {Pg  ILg)  ratio  versus  frequency  for  western  China  (top)  and  western 
United  States  {bottom).  Circles  show  average  ratios  for  earthquakes  and  asterisks  show 
average  ratios  for  nuclear  explosions.  Error  bars  are  for  1  standard  deviation.  Above 
about  3  Hz  separation  between  earthquakes  and  explosions  for  China  is  greater  than  for 
U.S.  and  scatter  appears  to  be  reduced. 
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